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Abstract 

A simple method for some class of inverse obstacle scattering problems is intro- 
duced. The observation data are given by a wave field measured on a known surface 
surrounding unknown obstacles over a finite time interval. The wave is generated 
by an initial data with compact support outside the surface. The method yields the 
distance from a given point outside the surface to obstacles and thus more than the 
convex hull. 
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1 Introduction and statements of the results 

The aim of this paper is to introduce a simple method for some class of inverse obstacle 
scattering problems in which some dynamical data over a finite time interval are used as 
the observation data. 

In order to explain the essence of the idea we consider two inverse obstacle scattering 
problems whose governing equations are given by the wave equations in three dimensions. 

1.1 Sound hard obstacles 

Let < T < oo. Let D C R 3 be a bounded open set with smooth boundary such that 
R 3 \ D is connected. Given / 6 L 2 (R 3 ) with compact support satisfying supp / n D — 
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let u = u(x,t) satisfy the initial boundary value problem: 

d 2 t u - Am = in (R 3 \D)x ]0, T[, 

^ = on&Dx]0, T[, 

dv (1.1) 
u(x,0) = inR 3 \ D, 

d t u(x,0) = f(x) inR 3 \D. 

Here we denote the unit outward normal to 3D by the symbol v. 

Let Q be a bounded domain with smooth boundary such that D C Q and R 3 \ f2 is 
connected. We denote the unit outward normal to dfl by v again. The dQ is considered 
as the location of the receivers of the acoustic wave produced by an emitter located at 
the support of /. In this paper first we consider the following problem. 

Inverse Problem I. Assume that D is unknown. Extract information about the location 
and shape of D from u on <9f2x]0, T[ for some fixed known f satisfying supp/ Pi O = 
and T < oo. 

This is a quite natural problem, however, to my best knowledge, it seems that no 
attempt has been done. Clearly the main obstruction is the finiteness of T and / is fixed. 
Note that u in (R 3 \ Q) x]0, T[ can be computed from m on 90x ]0, T\ by the formula 

u = z in(R 3 \n)x]0, T[ (1.2) 

where z solves the initial boundary value problem in R 3 \ Q: 

d 2 t z - Az = in (R 3 \ H) x ]0, T[, 

z — u on^fix ]0, T[, 

(1.3) 

z(x, 0) = inR 3 \fi, 

d t z(x,0) = f(x) inR 3 \H. 
Thus the problem can be reformulated as 

Inverse Problem I'. Extract information about the location and shape of D from u in 
(R 3 \ Q) x]0, T[ for some known / satisfying supp / fl Q = and T < oo. 

Now we state the result. Let B be an open ball with B fl Q = 0. Choose the initial 
data / G L 2 (R 3 ) in such a way that: 

(11) f(x) = a.e. x E R 3 \B, 

(12) there exists a positive constant C such that f(x) > C a.e. x G B or —f(x) > C a.e. 
x G B. 
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Let r > and v G i/ 1 (R 3 ) be the weak solution of 

(A - t 2 )v + f(x) = inR 3 . (1.4) 

This means that v satisfies 

- X7v ■ Vipdx - t 2 J vipdx = - J fipdx, Vy? G H^R 3 ). (1.5) 

The v is unique and is given by the explicit form 

v(x; t) = — j rf(y)dy, x G R 3 . (1.6) 

47T J B \x — y\ 

In this paper given two sets A and B we denote by dist (A, B) the distance between 
A and B: 

dist (A, B) = inf{|x-y| \ x e A,y e B}. 

If A or B consists of a single point, say B = {p}, we write dist (A, B) = c?a(p)- 
Set 

rT 

w(x;t)= e' Tt u(x,t)dt, x e R 3 \ fi, r > 0. (1.7) 

JO 

Our result is the following extraction formula from w and dwjdv on dQx ]0 T[ which 
can be computed from the data u in (R 3 \ Q)x ]0, T[. 

Theorem 1.1. If the observation time T satisfies 

T > 2dist(D,B) - dist (SI, B), (1.8) 

then there exists a r > such that, for all r > t 

fdv dw \ 
an \ou ou 



and the formula 



1 r ( dv dw \ 
lim — log/ \—w-—v\dS=-dist(D ) B) ) (1.9) 



Since dist (D, S) + ^/jftB|/47r coincides with the distance from the center of B to D, (1.9) 
yields the information about d£>(p) for a given point p in R 3 \J7. Therefore one can extract 
more than the convex hull of D. Note that we do not assume the special form of / except 
for the conditions (II) and (12). 

The restriction (1.8) on the observation time T is reasonable. Define the quantity 

l(dB,dD,dQ) =M{\x-y\ + \y - z\\x G OB ,y G 3D, zedVl}. 

This is the minimum length of the broken paths that start at x G dB and reflect at 
y G dD and return to z G d£l. We have 
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Proposition 1.1. 

2dist(D,B) - dist (SI, B) > l(dB,dD,dfl). 



Proof. One can find Xq G dB and yo £ ^.D such that | rc — 2/o I — dist (D, B). Let l(xo, yo) = 
{tx + (1 - t)y | < t < 1}. We see that l(x , y ) n <9fi 7^ 0. Let z G Z(x , yo) H <9fi. We 
have |xo — 20 1 > dist (Q, B). Thus 2dist (£>, _B) — dist (fl, 5) > 2|rr — yo| — |#o — zq\- Since 
l^o — ^0 1 — l^o — yo\ — 1 2/o — ^0 1 we have 2dist (D, £>) — dist (fl, B) > \xq — yo\ + |yo — za\- 
□ 

Therefore (1.8) ensures that T > l(dB,dD,dQ). This means that T is greater than 
the first arrival time of a signal with the unit propagation speed that starts at a point 
on dB at t — 0, reflects at a point on d£> and goes to a point on dQ. However, curiously 
enough in the proof of Theorem 1.1 we never make use of the finite propagation property 
of the signal governed by the wave equation. 

The procedure of extracting information about the location of D is extremely simple 
and summarized as follows. 

(i) Give an open ball B with B fl Q = 0. Using the initial data / satisfying (II) and (12) 
generate the wave field u. 

(ii) Choose a large T, say, such that T > 2 sup {\y — x\ | y G f2, x G B} — dist (fl, B) and 
measure u on dfl over the time interval ]0, T[. 

(hi) Compute the values of z in a neighbourhood of dfl relative to R 3 \ fl over ]0, T[ by 
solving (1.3). 

(iv) Choose a large r and compute u> and dw/du on <9f2 over the time interval ]0, T\ via 
(1.2) and (1.7). 

(v) Compute v and dvjdv on dfl via (1.6). 

(vi) Compute the quantity 



as an approximation of —dist (D, B). 

One choice of / gives one information about D by the procedure (i) to (vi). This 
means that we don't need to use many fs to get do(p) for a single p. This is the decisive 
character of our procedure. 

1.2 Penetrable obstacles 

Given / G L 2 (R 3 ) with compact support let u = u(x,t) satisfy the initial value problem: 

d 2 t u - V • 7VW = inR 3 x ]0, T[, 



where 7 = 7(2;) = (jij(x)) satisfies 

• for each i, j = 1, 2, 3 ^(x) = ^ji(x) G L°°(R 3 ); 

• there exists a positive constant C such that ^y(x)^ ■ £ > C|^| 2 for all £ G R 3 and a. e. 




w(x, 0) = inR 3 , 



(1.10) 



dtu(x,0) = f(x) inR 3 , 



x G R 3 . 
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This subsection is concerned with the extraction of information about discontinuity of 
7 from u on <9f2x]0, T[ for some / for a fixed T < oo. However, we do not consider the 
completely general case. Instead we assume: 

• there exists a bounded open set D with a smooth boundary such that 7(2;) a.e. x G 
R 3 \ D coincides with the 3x3 identity matrix J 3 . 

Write h(x) = 7(2;) — I 3 a.e. x G D. Our second inverse problem is the following. 

Inverse Problem II. Assume that both D and h are unknown and that one of the 
following two conditions is satisfied: 

(Al) there exists a positive constant C such that —h(x)£ ■ £ > |£| 2 for all (6R 3 and a.e. 
x G D] 

(A2) there exists a positive constant C such that h(x)£ • £ > |£| 2 for all £ G R 3 and a.e. 
x & D. 

Extract information about the location and shape of D from u on <9f2x]0, T[ for some 
fixed known f satisfying supp / fl Q = and T < 00. 

Note that w in (R 3 \ Q) x]0, T[ can be computed from -u on <9f2x]0, T[ by the exactly 
same formula as (1.2) and thus the problem can be reformulated again as 

Inverse Problem II'. Extract information about the location and shape of D from u in 
(R 3 \ Q) x ]0, T[ for some known / satisfying supp / fl Vt = and T < 00. 

Now we state our second result. 

Theorem 1.2. Assume that 7 satisfies (Al) or (A2). Let f satisfy (II) and (12) and v 
be the weak solution of (1-4) • Let T satisfies (1.8) and w be given by (1.7) with solution 
of (1.10). If (Al) is satisfied, then there exists a r > such that, for all t > r 



I i 

Jon \ 



dv dw \ 

—w — — v \ ab > 

OV OV I 



and the formula 



I f ( dv dw \ 
lim — log / — w-— v\dS = -dist(D,B), 
^2r s Jdn\dis dv J y ' 



is valid; if (A2) is satisfied, then there exists a r > such that, for all r > r 



Idn ( 



dv dw \ 

—w — — v ) ab < 

dv dv 



and the formula 



lim — lot 



is valid. 



Isakov [7] considered an inverse problem for the equation dfu — V • 7VM = in 
Ox ] — 00, T[ with the zero initial data u — when t < and the lateral Neumann data 
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du/dv = h on dQx ] — oo, T[, where f2 is the half-space x 3 < 0; 7 takes the value 1 for 
x G Q\D and a positive constant k for x E D. The D is given by a Lipschitz continuous 
function d on R 2 with d < as D = {x\ x 3 < d(xi, d 2 )}. 

His problem is to recover dD by measuring u on Tx]0, T[ for an arbitrary fixed 
nonempty open set T C <9f2 and known /i. Choosing /i as the lateral Neumann data of 
a special solution of the wave equation in the half-space, he showed that if k < 1, then 
the data u on Tx ]0, T[ uniquely determines the part of D within the set of all points 
x = (xi,X2,x 3 ) with x 3 > — T/2 and (xi,x 2 ) G Y. The condition k < 1 corresponds to 
(Al). It is an open problem whether or not the same conclusion holds in the case k > 1 
which corresponds to (A2). See also [8] for this point. 

Rakesh [12] considered an inverse problem for the equation 7<9 2 w — V • (jVu) = in 
R 3 x] -oo,T[ with u(x,t) = (t - x 3 ) 2 + for t « 0. Here (s) 2 = s 2 if s > and (s) 2 = if 
s < 0. The 7 takes 1 outside a bounded domain D with smooth boundary and a positive 
constant 1) on D. Thus the governing equation has a same constant speed inside 
and outside D. The data in his problem is the values of u on dQx] — 00, T[, where Q is a 
bounded open set of R 3 with smooth boundary and satisfies D C f2. He showed that; if 
T > 6 diam +inf :rg £) x 3 and D is strictly convex, then the data dujdv on (9f2x ] — oo,T[ 
uniquely determine D itself. 

Isakov employs a contradiction argument and his method starts with the uniqueness 
of the continuation of the solution of the wave equation and derives an orthogonality 
relation that was deduced by denying the conclusion. 

Rakesh's argument is also a contradiction argument and makes use of the uniqueness 
of the continuation of the solution of the wave equation. However, the main point is an 
analysis of the wave front set of u. See also [13] for other results. 

Unlike them we do not make use of the continuation of a wave field nor propagation of 
singularities argument. The method can be considered as an application of the enclosure 
method which was originally introduced for elliptic equations in [4, 3]. Recently in [5] 
the author found its application to inverse initial boundary value problems in one-space 
dimensional case for the heat and wave equations. In [6] we extended this method to the 
heat equation in two and three-space dimensional cases. Therein the initial data is zero 
and a special heat flux depending on a large parameter is used. 

1.3 Further remarks and construction of the paper 

Finally we comment on some results in the context of the Lax-Phillips scattering theory. 
Lax-Phillips in [9] established a relation between the support function of an obstacle and 
the right end point of the support of the scattering kernel which is the observation data 
in their theory. Since the support function gives the signed distance from the origin of 
coordinates to the support plane of the obstacle, the result means that one can get an 
estimation of the convex hull from the data. Note that the scattering kernel is written 
by using the scattered wave over the infinite time interval that is produced by a singular 
plane wave at t « far a way from the obstacle, and thus the data is completely different 
from ours. 

Majda[10] considered the singularity of the scattering kernel and clarified a relation 
between the support function and the right end point of the singular support of the back 
scattering kernel. In [11] a similar result for an obstacle with a finite refractive index is 
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given. The governing equation has the form a(x)d^u — Au = and a has a discontinuity 
across the boundary of the obstacle and takes 1 outside the obstacle. For other results 
including the Maxwell equations, hyperbolic systems, etc. we refer the reader to [10, 11] 
and references therein. 

A brief outline of this paper is as follows. Theorems 1.1 and 1.2 are proved in Subsec- 
tions 2.2 and 3.2, respectively. The key point of the proofs is to derive a lower estimate 
of the integral 

/.(£•-£•)-* ' ui > 

where v is the weak solution of (1.4). To establish the estimate we require some integral 
identities; these identities are found in Subsections 2.1 and 3.1. Using the identities, we 
show that, if T satisfies (1.8), then the dominant part in the lower estimate of (1.11) in 
Theorem 1.1 is essentially given by the integral of the square of v over D. We show that 
this last integral is comparable with e - 2r dist ( D, - B ) ignoring a multiplication of a power of 
r. This is stated in Subsection 2.2 and proved in Subsection 4.1. Note that in the proof 
of Theorem 1.2 instead of v the integral of |Vi> | 2 over D plays the same role and the 
corresponding estimate is stated in Subsection 3.2 and proved in Subsection 4.2. In the 
final section we give a conclusion of this paper and comments on further problems. 

2 The enclosure method for sound hard obstacles 

First we specify what we mean by the solution of (1.1). We follow the notion of the weak 
solution described on pp. 552-566 in [1] and use the notation therein. 

By Theorem 1 on p. 558 in [1], given u° G iJ x (R 3 \ D) and u 1 G L 2 (R 3 \ D) we know 
that there exists a unique u satisfying 

u G L 2 (0, T; H\R 3 \D)), u' G L 2 (0, T; H l (R 3 \D)), u" G L 2 (0, T; (H 1 (R 3 \ 13))') 

such that, for all G H\R 3 \ D) 

_ V«(x, t) ■ V(p(x)dx = a.e. t G ]0, T\ 

and u(x, 0) = u , u'(x, 0) = u 1 . In this section we say that this u for u° = and u 1 = f 
is the solution of (1.1). 

2.1 A basic identity 

Let u be the solution of (1.1). Define 

w(x;t)= / e~ Tt u(x,t)dt, xER 3 \D. 
Jo 

This w belongs to H 1 ^ 3 \ U). 



7 



Using integration by parts formula (Proposition 2 on p. 558 in [DL]), we see that, for 
all G if 1 (R 3 \ D) the w satisfies the equation 



/ _Vw ■ WSdx + / _(t 2 w - f)6dx 

JR3\D J-R?\D 

= -e~ TT [ _(u'(x, T) + tu(x, T)Udx. 



(2.1) 



/R 3 \D 

This means that, in a weak sense w satisfies 

(A - t 2 )w + f(x) = e- rT (u'(x, T) + tu(x, T)) in R 3 \ D, 

^ = 0ondD. 
av 

From_(2.1) for v G C™(R 3 \D), we have:(A - t 2 )w + f(x) = e~ Tl \u\x, T) + tu(x, T)) 
in R 3 \ D in the sense of distribution and hence Aw G L 2 (R 3 \ D). This yields that 
w G #f oc (R 3 \ D) and (A - t 2 )w + f(x) = e- TT (u'(x, T) + ru(x, T)) a.e. x G R 3 \ D. 

Now define dw/dv\dft as Vw\qq, ■ v G H l l 2 (dQ), where Vw\qq is the trace of Vw onto 

on. 

Let f G if 1 (R 3 ) be the weak solution of (1.4). For this v by the same reason as above 
we have v G ifj^R 3 ) and thus dv/dv\ 9n = Vv\ 9n ■ v G H l/2 (dVl). 
In this subsection we derive the following identity. 

Proposition 2.1. Let v G // 1 (R 3 ) be the weak solution of (1-4)- It holds that 



Idn ( 



dv dw \ 
av av I 



= \Vv\ 2 dx + T 2 \v\ 2 dx+ [ \V(w -v)\ 2 dx + T 2 I \w - v\ 2 dx 
Jd Jd Jr;^\d Jb?\d 

(2.2) 

+e~ rT f _(w - v)(u'(x, T) + tu(x, T))dx - e~ rT ( _(u'(x, T) + tu(x, T))vdx 

JR3\D JQ\D 

~~ / _f(w — v)dx— / fvdx. 

Jtt\D JD 

Proof. First we prove that 

Jon \av av J 

(2.3) 

= / — w — / _ f (w — v)dx — e~ rT \ _(u'(x,T) + Tu(x,T))v(x)dx. 
JdD av Jci\D Jn\D 

Let ip G H 1 ^ \ D) satisfy (p = on <9f2 in the sense of the trace. Since the zero extension 
of this ip belongs to if^R 3 \ D), it follows from (2.1) that 

/ _ Vw ■ V<p + I _{t 2 w - f + e- TT (u'(x, T) + ru(x, T))}<pdx = 0. (2.4) 
Jq\d Jn\D 



8 



Choose x £ C£°(R 3 ) such that x( x ) = 1 in a neighbourhood of Oil and x( x ) = in 
a neighbourhood of D. Since (1 — x) v \q\d vanishes in a neighbourhood of Oil, it follows 
that (2.4) is valid for ip — (1 — X^lnvD- 

On the other hand, since x v vanishes in a neighbourhood of D and w G ifj 2 QC (R 3 \ -D), 
integration by parts yields 



J _Vw ■ V{xv)dx 

/ — xvdS - / _(Aw)xvdx 
Jan Ov Jn\D 



[ ——vdS — [ _{t 2 w — f + e tT (u'(x,T) + ru(x,T))}xvdx. 
Jan Ov Ju\d 



! ^vdS = [ Vw-Vvdx + [ {r 2 w - f + e~ rT (u'(x,T) + Tu(x,T))}vdx. (2.5) 
Jan Ov Jn\D Jq.\d 



Ow 

"(IS - 

ln\D 

From this and (2.4) for <p — (1 — x) v \n\D we obtain 
Ow 

vdS = 

ln\D Jn\D ' 



Choose w G if x (f2) such that w = w in Q \ D. Note that w = w on Oil and OD in the 
sense of the trace. Since v G H 2 (il), we have 

Ov /" <9t> _ 



Jan Ov Jan Ov 

- / Avwdx + / Vv • Vwrfx 



- / _(r 2 f — f)wdx + / _ Vf • Vwdx 
Jn\D Jn\D 

+ / (r 2 t> — f)wdx + / Vi> • Vwdx. 
Jd Jd 



On the other hand we have 

Ov 



JdD 

that is 



/ Vi> • Vwdx = [ tt-w — ( (Av)wdx 

JD JdD OV JD 

= f ^—wdS — f (t 2 v — f)wdx, 

JdD OV JD 

f ^-wdS = f (t 2 v — f)wdx + / Vi> • Vwdx. 
JdD Ov Jd jd 



Therefore we obtain 

Ov r f f Ov 

/ _( r2-u — f)wdx + / 
ln\D Jn\D 

A combination of (2.5) and (2.6) gives (2.3). 



/ ^wdS = f _(r 2 v-f)wdx+ [ _Vv-Vwdx+ [ ^wdS. (2.6) 
Jan Ov Jn\D Jq\d JdD Ov 
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Write 

r dv f dv f dv 

/ —wdS = / %L( w -v)dS+ / ^vdS. (2.7) 

JdD OV JdD OV JdD OV 



i do dv JdD dv JdD dv 

It follows from (1.5) and the trace theorem that, for all <fi G i/ 1 (R 3 \ D) 

[ ^-(f)dS+ [ _ Vv ■ V(j)dx + r 2 / _v<f)dx = [ _f<f> 

JdD OV JH 3 \D JR3\D JR3\_D 

Combining this with (2.1), we obtain 
dv 

4>d$ — 

/R3\D v ' Jb?\D 



f -7r-4>dS — [ V (w — v) ■ V d>dx — r 2 [ (w — v)(bdx 
JdD dv Jh»\d Jn*\D 

/ (u'(x, T) + ru(x, T))6dx. 

JR3\D 



= e~* r 



/R3\D 

This means that w — v satisfies, in a weak sense 

(A-t 2 )(w-v) = e- rT (u'(x,T) + tu(x,T)) in R 3 \ D, 

9 l \ ^ an 

— — [w — v) — — — — onoD. 
dv dv 

Substituting w — v for <p in (2.8), we obtain 
dv 

[w-v)<M = 

/rs\d 1 v " Jr 3 \d 



I -^-(w-v)dS = I \W(w -v)\ 2 dx + T 2 [ \w-v\ 2 dx 
JdD dv Jr3\D Jr*\d 

+e~ rT ( _(u'(x, T) + tu(x, T))(w - v). 

JR 3 \D 



(2.8) 



(2.9) 



3 \ 

Now from this together with (2.3), (2.7) and the identity 

/ ^ v dS = r 2 [ \v\ 2 dx+ ( \Vv\ 2 dx- [ fvdx, (2.10) 
JdD dv Jd Jd Jd 

we obtain (2.2). 
□ 

In particular, choose / in such a way that supp/ Pi O = 0. Then (2.2) becomes 

Jan \dv dv J 

= \Vv\ 2 dx + r 2 \v\ 2 dx+ I \V(w -v)\ 2 dx + r 2 I \w - v\ 
Jd Jd Jr 3 \d Jr»\d 



'dx 



+e tT / _(w — v )(u'(x, T) + tu(x, T))dx — e rT / _{u'(x, T) + tu(x, T))vdx. 

(2.H) 

This is the basic identity for the sound-hard obstacles. 
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2.2 Proof of Theorem 1.1. 

Using the identity 

t 2 \w - v\ 2 + e- Tt {w - v)(u'(x, T) + tu(x, T)) 



-tT 



t(w -v) + (u'(x, T) + tu{x, T)) 

2r 



2 e -2rT 



At 2 



u'(x,T) + tu(x,T)\ 2 , 



we have from (2.11) 



/ \ tt-w — ^-v ) dS > t 2 f \v \ 2 dx 
Jan \av av J Jd 



-2tT 



/ _ k'fa;, T) + rit (:c, T)| 2 <ir — e rT / T) + tu(x, T))vdx. 

JR. 3 I) ./'' i> 



At 2 Jr3\d 
From (1.6) we have, for all x G R 3 \ B 



\v{x)\ < 



-Td B (x) 



li(b), |Vi;(x)| < 



-Tds(x) 



T + 



£i(B)- 



This gives ||u|| £ a (n \S) = C(e- rdist ( n > B )) and note that 

/ _\u'(x,T) + Tu(x,T)\ 2 dx = 0(t 2 ). 

JH?\D 

Using these estimates, we obtain 

f (?l w -^ v ) d s>T 2 [ ^dx + Oie-^ + Oire-^e-^^). 
Jan \ov av J Jd 

Here we state a key lemma whose proof is given in Section 4. 
Lemma 2.1. It holds that 

liminfT 6 e 2 ^( D < B ) / \v\ 2 dx > 0. 
r — >oo Jd 

Multiplying the both side of (2.14) by r 4 e 2rdlst (D ' B \ we have 

T 4 e 2.dist(D )S ) f (^ w _^ v ) dS > T S e 2Tdist(D,B ) f ^ 

Jan \av av J Jd 

+0(T 4 e" 2T(T " dist (D ' B)) ) + o(T 5 e~ r(T ~ 2dist ( D ' B )+ dist (^- B ))) 
and thus from (2.15) one gets 

Uminf r 4 e 2.dist (D,B) f ^ > Q 

r — >0 ° Jan \ov av J 
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provided if T > 2dist (D, B) - dist (ft, B). 

On the other hand, using (2.3), (2.7) and (2.10) one has 



Ian ( 



dv dw \ f f 

— w — —v\dS— / I Vf \ 2 dx + r 2 / M 2 cfa 

ov dv J Jd Jd 

dv 



(2.17) 



+ f ^ {w ^ v)dS + 0(re-^e- Tdist ^ B y). 

JdD dv 

From (2.13) we have, as r — > oo 

/ \Vv\ 2 dx + r 2 f \v\ 2 dx = 0(r 2 e- 2T(iist ^). 
Jd Jd 



(2.18) 



Concerning with the bound on the third term of the right hand side of (2.17), we have 
the following lemma. 



Lemma 2.2. It holds that, as r — ► oo 

dv 



dv 



Q^ T 2 e -rdlSt(D,B)^ 



and 



H- 1 / 2 (dD) 

,,,,,77, 0(rt ' ■ T 2 < ^ist(D,B)y 



li?!(R3\D) 

Proof. Set e — w — v. It follows from (2.9) that 



(2.19) 
(2.20) 



/ \Ve\ 2 dx + T 2 [ \e\ 2 dx = -e- rT [ (u'(x,T) + Tu(x,T))edx - [ ^-edS 
Jri\d Jb?\d Jk3\d JdD dv 



lR 3 \D 

and thus this yields 



ln 3 \D 



n 3 \D 



|V7,||2 _ , 2|| ||2 

I Ve llL 2 (R3\ J D) ~T T ll e llL 2 (R 3 \D) 



< e tT \\u'( • ,T) + tu{ ■ ,T) || L2 ( R 3\S)||e|| La(R 3 XS) + 



r dv 

JdD dv 



edS 



(2.21) 



The boundedness of the trace operator H l (R 3 \ D) — >■ H l / 2 (dD) yields 



dv 



edS 



8D dv 
< C 



< 



dv 



dv 



dv 



H- 1 / 2 (dD) 



\ e \dD\\H 1 / 2 (dD) 



(2.22) 



dv 



H- 1 / 2 (dD) 



lif!(R 3 \D)' 



where C is a positive constant independent of r. Now from (2.21) and (2.22) we obtain, 
for all 5 > 



:i-C5 2 )||Ve| 



L 2 (R 3 \D) 



^ + r 2 - C5 2 - 



Il 2 (r3d) 



5~ 2 
2 



< —e- 2rT \\u'( • ,T) + tu( ■ , T)|| 2 2(R 3 W) + C5- 2 



dv 



dv 



(2.23) 



ff-V2(aD) 
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Since v satisfies (A — t 2 )v = in D, for all \& G ^(D) we have 

^-mdS = r 2 [ v^dx + [ Vv ■ V^dx. 

dD OV JD JD 

The trace operator H 1 (D) — > H l l 2 (dD) has a bounded right inverse. This together with 
the identity above yields 



dv 



dv 



H- 1 / 2 (dD) 



<C(||V*;|| L 2 (d) +t 2 |M| l2(d) ), 



where C is a positive constant independent of r. It follows from (2.13) that this right hand 
side has the bound 0(r 

2 e -rdist(D,B)y Thug ( 2 19 ) ig vaM Now ^ 2 20 ) 

is a consequence 

of (2.19) and (2.23). 
□ 

Now we continue the proof of Theorem 1.1. 
It follows from (2.19), (2.20) and (2.22) that 

f ?±( w - v )dS = 0(rV^ T+dist ^ + r V 2 ^ dist < D > fl >). 

JdD OV 

From this, (2.17), (2.18) we obtain 

g 2rdist (D,B) f f^ w _^ v ) dS 

Jan \ov ov J 
= 0(r 3 e~ r(T ~ dist {D ' B)) + r 4 + re - r ( T - 2dist (A-B)+dist (n,B)y 
A combination of this and the estimate dist (D, B) > dist (f2, B) gives 



limsupr- 4 e 2Kiist ( D < B ) 



/ dv dw \ 
an \dv ov 



< oo 



(2.24) 



provided T > 2 dist (D,B) — dist (£l,B). Now the formula (1.9) is a direct consequence 

of (2.16) and (2.24). 

□ 



3 The enclosure method for penetrable obstacles 

First we specify what we mean by the solution of (1.10). By Theorem 1 on p. 558 in [1], 
given m° G if 1 (R 3 ) and u 1 G L 2 (R 3 ) we know that there exists a unique u satisfying 

u G L 2 (0, T; H^B 3 )), v! G L 2 (0, T; H\B?)), u" G L 2 (0, T; (^(R 3 ))') 

such that, for all <p G H l (R?) 

< u"{t), (f> > + [ 7(i)V«(x, t) ■ W<p(x)dx = a.e. t G ]0, T[ 

and u(x, 0) = u°, u'(x, 0) = u 1 . In this section we say that this u for u° = and u 1 = f 
is the solution of (1.10). 
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3.1 A basic identity 

Let u be the solution of (1.10). Define 

fT 



f 

w(x;r) = / e~ Tt u(x,t)dt, x E R 3 
Jo 



This w belongs to i/^R 3 ). 

From integration by parts (Proposition 2 on p. 558 in [1]) it follows that, for all <fi E 

J 7VW • V<pdx + J n3 (r 2 w - f)(f)dx = -e" rT J {u'(x, T) + ru{x, T))<pdx. (3.1) 

This means that in a weak sense w satisfies 

(V • 7 V - t 2 )w + f(x) = e- rT (u'(x, T) + tu(x, T)) in R 3 . 

By a similar reason as the sound-hard obstacle case we know that w E #f 0C (R 3 \ D). 
Since j(x) = 1 in H 3 \D, we define ^Vw-u\q^ as Vw\qq-u, where Vw\gn denotes the trace 
of Vw onto <9fl Note also that w satisfies (A — t 2 )w + f(x) = e~ TT (u'(x, T) + tu(x, T)) 
a.e. x E R 3 . 

In this subsection we derive an important identity. 
Proposition 3.1. Let v be the weak solution of (1-4) ■ It holds that 

/ {(Vf • v)w — {^Vw • v)v\dS 
Jan 

= — hVv-Vvdx + I 'iViw — v ) ■ V(w — v )dx + r 2 / \w — v\ 2 dx — / f(w — v)dx 
Jd Jr3 Jr3 Jn 

+e~ rT f (u'(x, T) + tu(x, T))(w - v)dx - e~ rT f (u'(x, T) + tu(x, T))vdx. 
jr» Jn 

(3.2) 

Proof. Using a similar argument for the proof of (2.5) and (2.6), we obtain 

/ (7VW • v)vdS = / 7VW • Vvdx + / {r 2 w — f + e~ TT (u'(x, T) + ru(x, T)}vdx 
Jen Jn Jn 

/ CVv ■ v)wdS = / Vf • Vwdx + / (t 2 v — f)wdx. 
Jan Jn Jn 



and 

From these we obtain 



/ {(Vf • v)w — {^Vw • v)v\dS 
Jon 

— / hVw ■ Vv dx — / f(w — v )dx — e~ rT \ (u'(x,T) + ru(x,T))vdx. 
Jd Jn Jn 



(3.3) 



Write 



- / hVw ■ Vvdx = - [ hVv ■ Vvdx - [ hV(w - v) ■ Vvdx. (3.4) 
Jd Jd Jd 
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(3.5) 



(3.6) 



Since v satisfies (1.5), from (3.1) we have, for all G iJ 1 (R 3 ) 

— / 7 V(u> — v) ■ Vddx — t 2 (w — v)d>dx 

= - / (J 3 - 7) Vv ■ V<pdx + e~ rT ( (u'(x, T) + tu(x, T))<f)dx. 

This means that the w — v satisfies, in a weak sense 

(V ■ 7 V - r 2 ) (w - v) = V • ( J 3 - 7) Vv + e- TT (u'(x, T) + ru(x, T)) in R 3 . 
Substituting w — v for in (3.5), we obtain 

/ (7 3 — 7 )Vf • V(w — v )dx = I {w — v) ■ V {w — v)dx 

iR 3 JR 3 

+r 2 / (w — v )(w — v )dx + e~ rT / T) + ™(x, T))(w — v)dx. 

J R 3 J R 3 

A combination of (3.4) and (3.6) gives 

— / hVw ■ Vvdx = — hVv ■ Vvdx 
Jd Jd 

+ / 7 V(u> — v) ■ V(iu — v)dx + t 2 (w — v)(w — v)dx 

JR 3 iR 3 

+e~ rT (u'(x, T) + tu{x, T))(w- v)dx. 

Now from this and (3.3) we obtain (3.2). 
□ 

In particular, choose / in such a way that supp / fl Q, — 0. Then (3.2) become 

/ {(Vi> ■ u)w — ( 7 Vw • u)v}dS 
Jan 

= — J hVv ■ Vvdx + y 3 7V (w — v) ■ V(w — v)dx + r 2 J ^\w — v \ 2 dx (3.7) 

+e~ TT [ (u'(x, T) + tu{x, T)){w - v)dx - e~ rT f (u'(x, T) + tu(x, T))vdx. 

This is our first basic identity which is useful in the proof of Theorem 1.2 under the 
assumption (A.l). 

Unfortunately, for (A2) this identity does not work. However, one can rewrite this 
by replacing the role of v and w in the proof of Proposition 3.1. More precisely, set 
/ = / — e~ rT (u'(x, T) +tu(x, T)). The points are: w satisfies V- 7 Vu> — r 2 w + f = inR 3 
and v satisfies V • I 3 Vv — r 2 v + f = —e~ rT (u'(x, T) + ru(x, T)) inR 3 . Thus changing the 
role of v and w in the proof of Proposition 3.1, we can easily obtain another expression 
of (3.2). 
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(3.8) 



Proposition 3.2. Let v be the weak solution of (1-4) ■ It holds that 

/ {(7VW • v)v — (Vf • v)w}dS 

Jdtt 

= I hVw-Vwdx + \ V(f — w) ■ V(f — w)dx + t 2 / \v — w\ 2 dx — / f(v — w)dx 
Jd jr 3 Jr. 3 Jn 

-e' TT f (u'(x, T) + tu(x, T))(y - w)dx + e~ rT f (u'(x, T) + tu(x, T))vdx. 

J R 3 J f2 

In particular, if supp / fl = 0, then (3.8) gives 

/ {(7VU! • z/)ti — (Vi> • v)w\dS 
Jan 

= J hVw-Vwdx + J V(v-w).V(v-w)dx + r 2 f 3 \v-w\ 2 dx (3.9) 

-e' rT [ (u'(x, T) + tu(x, T))(v - w)dx + e~ rT f (u'(x, T) + tu(x, T))vdx. 
jr. 3 Jn 

3.2 Proof of Theorem 1.2. 

First we consider the case when (Al) is satisfied. Using (Al) and the identity 

r 2 \w — v\ 2 + e~~ Tt {w — v)(u'(x, T) + ru(x, T)) 



o-tT 



t(w — v) H — - — (u'(x, T) + tu(x, T)) 
2r 



2 e -2rT 



4r 2 



u'(x,T) + tu(x,T)\ 2 , 



we have from (3.7) 



-2tT 



At 2 

Since we have 



/ {(Vv ■ v)w - hVw ■ u)v}dS > C [ \Vv\ 2 dx 
Jan Jd 

J \u'(x,T) + ru(x,T)\ 2 dx — e~ rT J (u'(x,T) + ru(x,T))vdx. 



(3.10) 



/ \u'(x,T) + ru{x,T)\ 2 dx = 0{t 2 ) 

JR 3 

and the estimate ||v||z,2(fi) = 0(e _rdlst (n ' B) ), it follows from (3.10) that 

= -rT^-rdist (n,B)y 



{(Vv ■ v)w - (7VW • u)v}dS > C / \Vvrdx + 0(e' 2rl ) + 0(re~ Tl e 
dn Jd 



Here we state a key lemma whose proof is given in the next section. 
Lemma 3.1. It holds that 



liminfT 4 e 2T ^ D ' B ) [ \Vv\ 2 dx > 0. 

T >oo J D 



(3.11) 



(3.12) 
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(3.14) 



Multiplying the both side of (3.11) by r ^e 2rdlst (D ' B \ we have 

r 4 e 2rdist(D, S) f { ^ v . u ) w ^^ w . u ) v}dS >r 4 e 2Tdist ^ [ \Vv\ 2 dx 

JdQ J D 

+0(r 4 e~ 2r(T ~ dist ( D - B ))) + 0(r 5 e~ r(T ~ 2dist ( D > s )+ dist ( n > B ))) 
and thus from (3.12) one gets 

liminf r 4 e 2rdist(D ' B) f {(Vv ■ v)w - (7W? ■ v)v}dS > (3.13) 

provided if T > 2dist (D, B) - dist (fl, B). 

On the other hand, using (3.3) and (3.4), one gets 

J {(Vf • v)w — {jWw • v)v}dS 

= - I hVv Vvdx - [ hV{w-v)- Vvdx + 0(re- rT e- rdist ^). 
Jd Jd 

From (2.13) we have, as r — ► oo 

/ \Vv\ 2 dx = 0{r 2 e' 2T ^ t{ - D ^). (3.15) 

JD 

Concerning with the bound on the second term of the right hand side of (3.12), we have 
the following lemma. 

Lemma 3.2. It holds that, as r — > oo 

\\w - u||hi ( r3) = 0(re- rT + T e^ dlst ^) 

Proof. Set e — w — v. From (3.6) we have 

- J hVv ■ Vedx - J 7Ve • Wedx - r 2 J \e\ 2 dx = e~ rT f R3 ( u> ( x , T ) + tu(x, T))edx 
This yields 

C||Ve|| 2 2( R 3) + 7" 2 |HIl2(r3) 

< e~ rT ||w'( • ,T) + tu( ■ ,T)|| L 2( R 3)||e|| L 2 (R 3) + j^hVv - Vedx 
and thus one gets, for all 5 > 

(r 52 Wn 2 + 52 W 

I c ~ I H Ve llL2(R3) + I r - — I ||e|| L 2 (R 3) 



(3.16) 



r)~ 2 

< : d —e-^ T \\u" 
~ 2 



'(•,T) + ™(-,T)|| L2(R3) + —\\h\\ L ^ (D) \\Vv\\ L 2 [D) . 
Now a combination of this and (3.15) yields (3.16). 
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□ 

A combination of (3.15) and (3.16) gives 



Jd 



j hV(w -v) -Vvdx = 0(T 2 e~ T 



(T+dist (D,B)) , 2 -2rdist (D,B) 




/ ((Vf • v)w — (7VW • v)v) dS 
Jan 



= 0(r 2 e" r(T ~ dist (D ' B)) + r 2 + re - T ( T - 2dist (A-B)+dist (n,B)^ 
This together with the estimate dist (D, B) > dist (f2, B) yields 



limsupr- 2 e 2rdist(D ' s) [ ((Vv ■ v)w - (7VW • v)v) dS < 00 



r — >oo JdU 



(3.17) 



provided T > 2 dist (D, B) — dist (fl, B). Now the conclusion of Theorem 1.2 is a direct 

consequence of (3.13) and (3.17). 

□ 

Finally we give a comment on the case when (A2) is satisfied. In this case we make 
use of (3.9) instead of (3.7). A combination of the well known inequality (see [2]) 



(7(2;) - h)Vw ■ Vw + V(w - v) ■ V(w -v)> (7(0:) - h)^{x)~ 1/2 Vv ■ ^(x)~ 1/2 Vv 
and (A2) yields that there exists a positive constant C such that 



which corresponds to (3.10). Applying the arugument for the proof of (3.13) to this right 
hand side we obain 



T > °° JdVL 

Since (3.17) is valid also for case (A2) we obtain the desired conclusion. 
□ 

4 Proof of Lemmas 2.1 and 3.1 
4.1 Proof of Lemma 2.1. 

Choose points x e 3D and y G dB such that dist (D,B) = \x — y \. Since we have 
assumed that dD is smooth, one can find an open ball B' such that B' C D and x G 



h(x)Vw ■ Vw + V(w - v) ■ V{w -v)> C\Vv 



for a.e. x G D. This together with (3.9) gives the lower estimate 





f \u'{x,T) + ru{x,T)\ 2 dx + e~ TT f {u'{x,T) + ru{x,T))vdx 
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dB' fl dD. Since dist (B',B) = \x — yo\, it suffices to prove (2.15) in the case when 
D = B'. 
Write 



I \ 2 r g-^(|j/l-^l + l?/2-a;|) 

/bxb |x - y\\\x - 2/2 1 
It follows from the assumption on / that 



/ 1 \ A r g-TU»l-x|+|j/ 2 -x|; 

= j- / i n rf(yi)f(y2)dyidy 2 . 

\A-kJ Jbxb \x — y 1 \\x — yo\ 



v(x) 2 >C 2 I(x,r) 2 (4.1) 

where 

1 r e ~r\y-x\ 

H x , r ) = t~ / i r^' x G 

47T Jb \v — x\ 



We denote by p and 77 the center and radius of B, respectively. Using the polar 
coordinates centered at x, one can write 

B = {y = x + ru\ujE S(x, B), r + (u) < r < r^{uj)} 

where 



S(x, B) = {lo G S I to ■ (p — x) > y \x — p\ 2 — r] 2 } 

and 



r ± (u;) = uj ■ (p — x) =F \J iy> • (p — x )) 2 — \x — p\ 2 + rf. 
This together with r + {uj) > dist (B', B) yields 

rr-{u>) 



I(x, r) > dist (£', B) f duj e~ Tr dr 

JS(x,B) Jr+(oj) 



and thus we have 



Tl {x, r) > ^ig^gj I f e- Tr+ ^du - \ e-^^duj) . (4.2) 

47T \Js(x,B) Js(x,B) J 

Define d B {x) = inf{|rr — y\ \ y G B}. First we give an estimate for the second integral 
in the right hand side of (4.2). Since ds(x) = \x — p\ — rj and r~(oj) > \J\x — p\ 2 — rf, we 
have 

/ e~ Tr ~ {u)) duj < [ e -^V^)V\ x -p\+^duj 

JS(x,B) JS(x,B) 



IS(x,B) JS(x,B) 



Here note that 



2r]Jd B {x) 

^d B {x)^\x -p\+ri- d B (x) = == > 0. 

\x -p\ + i] + Jd B {x) 



Thus we obtain 

e TdBi - x) ( e- Tr ~^duj < 47re- Ar (4.3) 

JS(x,B) 
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where 

A = inf jj^/fgg 

XeB ' ^\x-p\+T] + y/d B {x) 

The points are: A is positive and independent of x G B' . 

Next we consider the first integral of the right hand side of (4.2). The surface S(x, B) 
has the parameterization: 

p — x 



ou(r, 9) — r cos 9a(x) + r sin 9 b(x) + , jVl — r 2 , (4.4) 

p — x\ 

where < r < i]/\x — p\ and < 9 < 2tc; a(x) and b(x) are unit vectors perpendicular 
each other and satisfy a(x) x b(x) — {p — x)/\p — x\. Since duo = (r/yl — r 2 )drd9, one 
can write 

Js(x,b) Jo Jo v 1 — r 2 

Here we note that r + (u(r, 9)) is independent of 9. In fact we have 



r + (co(r, 9)) — \p — x\ Vl — r 2 — \Jrf — r 2 \p — x\ 2 . 

Thus this yields 

/ e~ Tr+ ^duj > 2tt [ V/lP ~ Xl e' rr+ ^ r ^rdr. (4.5) 

Js{x,B) Jo 

Since d B (x) = \p — x\ — rj, we obtain 

r + (uj(r, 0)) — d B (x) — \p — x\y/\ — r 2 — \Jrf — r 2 \p — x\ 2 — \p — x\ + rj 



\p — x\(y/\ — r 2 — 1) + (rj — \jrj 2 — r 2 \p — x 



2\ 



r 2 \p — x\ 



(4.6) 



jj _|_ ^2 _ r 2|p _ x \2 1 + y/\ — 



r 2 \p — x\ r 2 \p — x\ 
< = ' < — — - 



fj + yr/ 2 — r 2 \p — x\ 2 V 
Now set M = sup(|p — x\/rj). Since r + (cu(r, 0)) — d B (x) > 0, from (4.6) we obtain 

x&B' 

m/\x-p\ rv/\p~A 



[ e- Tr+( ^rdr > e' rd ^ f ' r^^rdr 

Jo Jo 

e ~Td B {x) 



2tM 



.(1 _ e -rB( V /\p-x\)^ (4,7) 



-rd B (x) 

> (l- e - T / B ). 
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From this, (4.2), (4.3) and (4.5) we can conclude that: there exist r > and C > 
independent of x G B' such that, for all r > t 

tV^/^t) > C". 

This together with (4.1) yields 

r 4 / v(x) 2 dx > (CC'f ( e- 2TdB{x) dx = (CC') 2 e 2Tr > [ e~ 2r ^dx. 

J B' JB' J B' 

We have already known (Proposition 3.2 in [6]) that 

liminf rV rd B' (p) / e' 2rlx ~ pl dx > 0. (4.8) 

T >OC J B , 

This yields 

liminfr 6 e 2r( ^'(p)-") / v (x) 2 dx > 0. 

T >00 J B , 

Since ds'ijp) — V — dist (£>', £?) we conclude that (2.15) is valid. 
□ 

4.2 Proof of Lemma 3.1. 

We employ the same notation used in the proof of Lemma 2.1. It suffices to prove (3.12) 
in the case when D = B' . Since 

Vv{x) = -I j B 1 \l^~f e-^f{y){x - y)dy, 

we have 

(47r) 2 |Vt;(x)| 2 

, (4-9) 
= / K(x, yi )K(x,y 2 )e-^-y^-y^f( yi )f(y 2 )(x - Vl ) ■ (x - y 2 )d yi dy 2 , 

JBxB 

where 

v( \ l + r\x-y\ 2 
K(x,y) = 



\x-y\ 



3 



Define 

Since it holds that (x — y±) ■ (x — y 2 ) > for all y±,y 2 G B and x G R 3 \ B and 
f{yi)f{V2) > C 2 a.e. 2/1,3/2 G B, from (4.8) we obtain 

(47r) 2 |Vt;(x)| 2 

> C 2 f K(x, Vl )K(x, y 2 )e-^ x - y ^ x - y ^(x - Vl ) ■ (x - y 2 )dy x dy 2 

JBxB 

(4.10) 

> C 2 (l + rdist (D, B) 2 ) 2 ( f ~ Vl ]' {X - y A e -<\ x -^ + \ x -y^d yi dy 2 

Jbxb \x — yi\ 6 \x — y 2 \ 6 

= C 2 (l + rdist (D, B) 2 ) 2 \J(x, r)| 2 . 
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One can write 



J(x, t) — uodto I 

JS(x,B) Jr 



r (w) 

a; da; / e rr <ir 

S(x,B) Jr+(u) 



(4.11) 



-[/ e- Tr+{u}) uuduu- [ e~ Tr ~ {u}) ujduj 

T \Js(x,B) JS(x,B) j 



Since r + {uj) > ds^x), we have 



JS{x,B) 



< 4ire- TdB(x) . 



The second integral in the right hand side of (4.11) has the bound 0(e rdB ^e At ). These 
give 



r 2 \J(x,r) 



f e- Tr+(u}) ujduj 

JS(x,B) 



+ 0(e~ AT e~ 2TdB ^). 



(4.12) 



Using parameterization (4.4) of S(x,B), one has 



- Tr+(uj) udu 



S(x,B) 



•2vr rv/\x-p\ 

d9 
o Jo 



-Tr+(u>{r,0)) 



u(r,9) 



rdr 



Since r + (o;(r, 9)) is independent of 9 and 



r 2ir 

/ u(r,9)d9 
Jo 



p — x 
\p — x\ 



Vl -r 2 27r, 



we obtain 



Js(x,B) \p — X\ Jo 



lS(x,B) \p ■ 

From this together with (4.7) and (4.12) one can conclude that: there exist tq > and 
C > independent of x e B' such that, for all r > r 

r 4 e T2dB ^\J(x : r)\ 2 > C . 

Thus from (4.10) we obtain 

t 2 I \Vv(x)\ 2 dx> C" I e- 2TdB{x) dx = C"e 2Tr < [ e^^dx, 
Jb' Jb' Jb' 



where C" is a positive constant. Hereafter using (4.8), we obtain (3.12). 
□ 



5 Conclusion and further problems 

In this paper we introduce a simple method for some class of inverse obstacle scattering 
problems that employs the values of the wave field over a finite time interval on a known 
surface surrounding unknown obstacles as the observation data. The wave field is gener- 
ated by an initial data localized outside the surface and its form is not specified except 
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for the condition on the support. The method yields information about the location and 
shape of the obstacles more than the convex hull. 

• It would be interesting to apply the method presented in this paper to other time 
dependent problems in electromagnetism, linear elasticity, classical fluids etc.. Those 
applications belong to our future plan. 
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